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ABSTRACT 

Swift XRT observations of the Hi line source VIRGOHI 21 were performed on 
April 22 and April 26, 2008 for a total exposure time of 9.2 ks. This is the first 
pointed X-ray observation of VIRGOHI 21, a putative dark galaxy in the Virgo 
cluster, and no photons were detected from this source. The non-detection of 
extended X-ray emission within the angular extent of the Hi source corresponds 
to a 99% confidence upper limit of 2.1 x 10~ 14 ergs cm~ 2 s _1 in the 0.3-2.0 keV 
band. The equivalent upper limit to the amount of diffuse hot gas associated 
with VIRGOHI 21 is in the range 4 x 10 7 — 2 x 10 8 M for a hot gas temperature 
between 0.1 and 1 keV. The non-detection also corresponds to a 99%-confidence 
upper limit on the flux from a point-like source of 8 x 10 -15 ergs cm -2 s _1 in 
the 0.3-2.0 keV band. We discuss the constraints on the nature of VIRGOHI 21 
imposed by these observations and the theoretical implications of these results. 



Subject headings: dark matter 
galaxies 



galaxies: individual (VIRGOHI 21) — X-rays: 



Introduction 



Hi surveys of the Virgo Cluster flDavies et al.l 12004 12009 : iMinchin et all 120051 . 12007 



Haynes et al.l 120071 ) have revealed a remarkable Hi cloud, VIRGOHI 21, with apparent 
evidence of rotation yet unaccompanied by any optical emission. The Hi mass (Mhi ~ 
3 x 10 7 M ), inferred circular velocity (V ~ 100 km s _1 ) and physical size (R ~ 8 kpc) are 



department of Physics, University of Alabama, Huntsville, AL, USA 
2 NASA Marshall Space Flight Center, Huntsville, AL, USA 

3 Universities Space Research Association, NASA Marshall Space Flight Center, VP62, Huntsville, AL, 
USA 

4 Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, MS-67, Cambridge, MA 02138 



-2- 



all typical of luminous galaxies of dynamical masses 10 10 — 10 11 M yet deep optical follow-up 
observations provide an upper limit of onl y 31.1 / mag arcsec -2 , corresponding to a stellar 



content of 2xl0 5 L Q flMinchin et al.l 120071 ). Thus, VIRGOHI 21 is currently considered the 



prototypical dark galaxy candidate: a dynamically-massive dark-matter-dominated galaxy 
with a rotating gaseous disk but devoid of stars. 

If VIRGOHI 21 is a massive dark galaxy, then it has important implications for the 
theory of structure formation. Th e high ly-successful cosmological hierarchical clustering 
model originated by IWhite fc Reed (119781 ) is based on the idea of accretion and merging 
of cold dark matter (CDM) from small sizes (~10 8 M©) up to the large galaxy clusters 
observed at the present epoch. The CDM model predicts a large number of ^10 8 M Q dark 
matter (DM) halos sho uld remain unassimilated following hierarchical growth an d thus still 



be present today (e.g., iMoore et al.lll999l ; iKlypin et al.lll999l ; iDavies et al.l 120061 ) . To date 



however, the number of observed low-mass galaxies falls far short of these predictions. This 
is tru e locally, wh ere only a handful of satellites surround the massive Local Group spirals 



Mateo 



19981 ). and is manifest in the fai nt end of the luminosity function that s hows a 



(e.g. 

slope more shallow than ACDM predictions (iBlanton et al.ll200ll ; iBenson et al.l 120031 ). This 
shortcoming of the ACDM model can be overcome either by preventing small DM halos from 
forming in the first place or by suppressing star formation so that they remain optically 
faint. The latter possibility implies a population of low- mass DM halos in which the baryon 
content remains optically-dark but still visible chiefly through Hi gas emission. T he rarity 
of isolated extragalactic Hi clouds lacking optical counterparts (jDoyle et al. 2005 ) suggests 
that VIRGOHI 21, if truly massive, is a rare case in favor of this scenario. 



An alternative to the dark galaxy inte rpretation for VIRGOHI 21 (IBekki et al.l 12005 



Haynes et al.l 120071 : iDuc fc Bournaudl 120081 ) is that it is simply a part of an elongated Hi 
dominated tidal tail emanating from the nearby luminous galaxy NGC 4254 (M 99). The tidal 



tail w as caused either by a fly-by encounter with another (massive) galaxy (jVollmer et al. 



20051 ) or by ram pressure stripping within the Virgo environment. Recent observations 
show the Hi gas in the tail extends from NGC 4254, located ~120 kpc to the south o f 
VIRGOHI 21, to about 130 kpc beyond VIRGOHI 21 to the north iHavnes et al.l l2007h . 
The peculiar velocity structure in the vicinity of VIRGOHI 21 is, in this case, ascribed to 
projection effects on the streaming motions of the debris tail instead of to rotation within a 
gravitationally bound disk. In this scenario, the estimated Hi mass is close to the total mass 
of the object; there is no massive halo accompanying VIRGOHI 21. 

In § [2] we describe the Swift observations leading to upper limits of the X-ray emission 
from VIRGOHI 21, and in § [3] we provide our interpretation of the Swift non-detection. In 
this paper we assume a distance to the source of D = 16.5 Mpc corresponding to the mean 



- 3- 



of the Virgo Cluster (IMei et al.l 120071 ) and a scale of 4.9 kpc per arcminute. 



2. SWIFT observations of VIRGOHI 21 and data analysis 

The VIRGOHI 21 field was observed by Swift on 2008 April 22 and 26 with the XRT 
operated in photon counting mode. The XRT has a very low detector background, which 
makes it especially suitable for observations of faint diffuse X-ray sources. These data are 
the first pointed X-ray observations of VIRGOHI 21. 

The data were screened for bad pixels and other detector artifacts, and only photon 
grades 0-12 were included in the analysis. The screening resulted in 9.2 ks of clean data. In 
this analysis we considered photon energies in the 0.3-7 keV range, where the calibration of 
the XRT is better understood. The data reduction was performed using xselect and the 
FT00LS, and the spectral analysis with XSPEC. 

VIRGOHI 21 lies about 1 Mpc to the NW of the center of the Virgo Cluster at 
R.A.=12 h 17 m 53.6 s decl.=+14°45'25" (J2000.0). The a ngular size of VIRGOH I 21, 3.5'x 1.375' , 



lies w ell within the 24'x24' field of view of the XRT (jCusumano et al.ll2006l ; iRomano et al. 
20051 ). X-ray images of the VIRGOHI 21 field are shown in Figure dj in which the gray 



contour are from the Hi observations of iMinchin et al.l (120071 ). The upper part of the Hi 
feature was interpreted by lMinchin et al.l (120071 ) as a rotating disc, and the lower part as an 
Hi bridge connecting VIRGOHI 21 with NGC 4254, located to the south and outside the 
Swift field of view. 

As with normal galaxies, VIRGOHI 21 was expected to be a source of both diffuse 
and point-source X-ray emission. Even in the absence of stars, gas falling into a galaxy's 
gravitational potential from the intergalactic medium should be heated to roughly the (X- 
ray-emitting) virial temperature. In addition, accretion onto a central supermassive black 
hole from the gas reservoir in VIRGOHI 21, even at a fraction of the Bondi rate, should 
appear as a localized source of X-radiation though not necessarily optically bright. No X-ray 
emission associated with VIRGOHI 21 was detected. An unrelated weak (41 counts) point- 
like source was detected near R.A. = 12 /l 17 m 23 s , decl.=14°39'29" (Figured]). As no object is 
known within 1' of this location, we regard this as the discovery of the X-ray source SWIFT- 
XRT J121723+1439.5. In the following, we place upper limits to the presence of diffuse and 
point source emission from VIRGOHI 21. 



-4- 



2.1. Diffuse X-ray emission 

In order to place upper limits to the diffuse soft X-ray emission from VIRGOHI 21, we 
accumulated X-ray photons from a 3.5' x 1.375' region located at the position of the putative 
Hi disc (black box in Figured]). We detected a total of 8 counts in the 0.3-2 keV band, and 
no counts in the 2-7 keV band. We estimate the number of background counts, first by 
extracting the data from the entire field of view, excluding a region of enhanced emission 
to the south (Figure [I]) and regions close to the detector boundaries. The region included 
encompassed 337.3 square arcminutes. We found 623 counts in the 0.3-2 keV band and 319 
counts in the 2-7 keV band. Ignoring the possibility that some of these events might be 
vignetted by the telescope one therefore expects in the region of interest at least 8.9 ± 0.4 
counts in the soft band and 4.6 ± 0.3 in the hard band. We note that it is mildly interesting 
and somewhat improbable that no counts were measured in the hard band from the region 
containing VIRGOHI 21. 

We use the soft X-ray band data in order to set upper limits as to source detection, since 
this is the band in which the X-ray emission from diffuse gas (and point sources) is expected. 
Since the measured background predicts at least a mean of 8.9±0.4 counts and we detected 
8 counts it is therefore evident that the VIRGOHI 21 region is not an X-ray emitter. To 
determine upper limits we increased the expected number of counts in the soft band by 10% 
to 9.8 counts to account for the possibility that some of the measured background events arise 
from faint, unresolved X-ray sources. We determine 99%-confidence upper limits assuming 
18 counts as detection. 



2.1.1. Uniform distribution 



We assume that the Hi disk is seen edge-on, and thus the emitting volume is V = ttR 2 L, 
where R = 1.75' and L = 1.375'. We use an optically-thin plasma of uniform density and of 
Solar abundance (apec in XSPEC), a Galactic Hi col umn density of Nh 



2.7 x 10 20 cm" 2 



( IDickey fc Lockmanl Il990l ; iKalberla et al.l 120051 ) with iMorrison &: McCammonl (119831 ) cross- 
sections (wabs in XSPEC), and calculate the emission measure required to achieve 8.1 source 
counts as a function of plasma temperature using XSPEC. The emission measure calculated 
by XSPEC is given by the model normalization, K, as 



K 



10 



-14 



4tiD 2 a 



n e n,HdV, 



(1) 



in which n e and are, respectively, the electron and hydrogen number density and Da the 
angular size distance (for z<l). From this, we calculate the gas density, and thus the gas 



- 5 - 



mass. The results are shown in Figured showing that these 9.2 ks Swift XRT observations 
set upper limits of the order of 10 8 M Q . These upper limits correspond to a flux of 2.1 x 10~ 14 
ergs cm -2 s _1 in the 0.3-2 keV band. 



2.1.2. (3 model distribution 

A putative hot halo may be maintained in hydrostatic equilibrium by the gravitational 
potential of the galaxy. In this case, the gas density will decrease with radius. We therefore 
consider a more realistic model of the gas density given by the (3 model, 

n e = n eQ \l + —\ (2) 



as often observed in galaxy clusters (jCavaliere fc Fusco-Femianolll978l ). The parameters r, 



and (3 describe, respectively, the radius within which the density is approximately constant, 
and the decrease of the density at large radii. For this interpretation, we use a spherically 
symmetrical model in which the galaxy has a radius of R = 1.75'. In this region, the 
background in the 0.3-2 keV band has a mean of 16.3 counts, and we thus place the 99% 
confidence upper limit at 27 total counts (or 10.7 source counts), following the method 
described in § 12.1.11 Using the (3 model (Equation [2]) into Equation [Q, we obtain upper 
limits to the central gas density: 

( KWHkD* \ * 
n e0 = ~ A — TTTFTi \ ( 3 ) 



A-nrll^R/r,. 



in which the integration of the density profile has led to the integral I\(R/r c ) = f^ Tc x 2 (l + 
x 2 )~ 3l3 dx. We use two r eference values for the index {3: a value of /?=2/3, as usually measured 



in galaxy clusters (e.g.. lLaRoque et al.ll2006l ). and a shallower slope of /3=l/3, which may be 
typical of lower-mass halos (e.g., galaxy groups and galaxies). We find analytical solutions 
for the integral, I\{n) — n/(2 + 2n 2 ) + arctan(n)/2 for (3=2/3, and Ii(n) = n — arctan(n) 
for 0=1/3. 



We then estimate the upper limits to the gas mass by integration of the density profile, 
M gas = [i e m H 







L neo ( i_ 






r 2J 



-1/3 



= 47Tfi e m H n e0 r 3 I 2 {R/r c ), (4) 

in which the (3 model profile has led to the integral I^iR/r^) = f^ Tc x 2 {l + x 2 )~ 3/3 ^ 2 dx. This 
integral is ^2( n ) = n — arctan(n) for (3 = 2/3, and ^{n) = hi(x + yl + x 2 ) for (3 = 1/3. 
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We estimate upper limits to the gas mass as function of plasma temperature using 
Equations [3] and HI using a fiducial value of the core radius of 0.5' (Figure [2]). We determine 
that the upper limits to the gas mass depend weakly on the choice of the core radius, given 
that varying the fiducial value by ±100% results in variation of less than 10%. This model 
results in upper limits to the gas mass that are in the range of M gas ~ 4 x 10 7 — 2 x 10 8 
Mq. Upper limits to the gas mass are therefore similar to those obtained from the simple 
uniform density model of Section 12.1.11 

2.2. Point source X-ray emission 

We set upper limits to the luminosity of an unobserved point source anywhere within 
the 3.5' x 1.375' region of interest. The upper limit takes into account that there are 68 
statistically-independent SWIFT spatial resolution elements within this region, given that 
the half-power diameter of the XRT point spread function at 1.5 keV is 18". The measured 
background, per 18"-diameter spatial resolution element is 0.131, or 0.144 using a +10% 
allowance for vignetting. Thus, considering the Poisson distribution with a mean of 0.144 
and a total of 4 counts, there is a probability of only 1.6 x 10~ 5 of measuring 4 or more counts 
in just one of these resolution elements. Of course we would have been happy to measure 
4 or more counts in any one of the 68 elements so the joint probability is correspondingly 
larger. Thus, we may state, with approximately 99.9%-confidence that we have not detected 
a a source of 4 counts or brighter from this region. 

We consider a power-law emission model with photon index V = 2.0 and determine the 
upper limits to the source luminosity as function of the absorbing Hi column density using 
XSPEC. The upper limits are shown in Figure [21 and are in the L = 3 x 10 38 — 10 40 ergs s~ x 
range. This upper limit corresponds to a flux of 8xl0~ 15 ergs cm -2 s _1 in the 0.3-2 keV 
band. 



2.3. UV and IR emission 

For completeness, we also analyzed extant non- X-ray images of VIRGOHI 21. 

Swift UVOT images using the UVM2 filter (effective wavelength 2231 A) were obtained 
concurrently with the XRT data. No sources were detected within the angular extent of 
VIRGOHI 21. The FT00LS utility uvot source was used to estimate a limiting magnitude 
for the UV emission from the 3.5'x 1.375' region enclosing VIRGOHI 21. This limit does not 
improve our estimates of the upper limit to the mass of diffuse hot gas nor to the flux from 
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a point source as estimated above from the X-ray data. 

We also analyzed archival Spitzer data (Program ID 30725) which included all four 
IRAC and three MIPS spectral bands. Again, no extended emission was detected within the 
Hi extent of VIRGOHI 21. We estimate an upper limit to the mass of dust associated with 
VIRGOHI 21 to be 10 4 -10 6 M depending on the dust temperature assumed. 



Discussion and conclusions 



A search for X-ray emission from the Hi line source VIRGOHI 21 using the Swift XRT 
obtained only upper limits. The non-detection indicates that any diffuse X-ray emitting gas 
at kT = 0.1 — 1 keV must be less massive than 4 x 10 7 — 2 x 10 s M©, depending on the 
gas temperature and its spatial distribution. Given that the Swift XRT detector has no 
effective area approximately below 0.1 keV, it is possible that significant amounts of gas at 
sub-virial temperature have gone undetected by this X-ray observation. Similarly, any X-ray 
point source associated with VIRGOHI 21 will have a luminosity lower than 5 x 10 38 — 10 40 
ergs s _1 , depending on the absorbing column density of gas along the line of sight. 

How do these limits compare with expectations? There are well-established correlations 
between optical and X-ray light from galaxies: Typical fx/ fo ratios for normal galaxies 
range from 10~ 4 to 10~ 2 and for galaxies hosting AGNs the range extends up to unity (e.g., 
Kim et al.ll2006l ). For the dark galaxy VIRGOHI 21, these ratios predict an X-ray flux of 
10~ 17 to 10~ 15 for a normal galaxy up to 10~ 13 if an AGN is present. Of course, these fx/ fo 
correlations reflect the causal link between stellar processes and X-ray light. They are not 
expected to hold in the absence of stars. 

Instead, our expectation that V IRGOHI 21 should be a source of diffuse X-ray emission 
stems from the fundamental picture (IRees fc Ostrikerlll977l ) that gas falling into a DM halo is 
first heated to the virial temperature. The accumulation o f hot gas over a Hubbl e time could 
be substantial for a massive DM halo like VIRGOHI 21 Jwhite fc FrenkHl99l[ ). We would 
expect a gas virial tem perature of approxim ately 0.4 keV based on the observed kinematics 
of VIRGOHI 21 (from iMinchin et al.l 120071 ) and thus our best guess to the amount of hot 
gas associated with VIRG OHI 21 is < 5 x 10 7 Scaling from simulations of gas infall 
onto stable disk galaxies (IBenson et al.ll2000l ; iToft et al.l 120021 ). the X-ray luminosity of this 
gas should be up to 5xl0 38 ergs s _1 for rotation speeds similar to that of VIRGOHI 21. 

The above estimates are based on the assumption of an isolated 'quiescent' massive 
galaxy accreting slowly (i.e., with a long cooling time) from the surrounding intergalactic 
medium. In the particular case of VIRGOHI 21, in contrast, a considerable gas reservoir 
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exists within the tidal tail emanating from NGC 4254 that should lead to a relatively high 
gas accretion rate and hence high X-ray luminosity. Furthermore, turbulence and shocks 
following the interaction between VIRGOHI 21 and NGC 4254 would very likely heat some 
portion of the disk Hi gas component of VIRGOHI 21 to X-ray emitting temperatures. We 
therefore consider the above estimates to be minimal expectations. Of course, any prediction 
of the amount of hot gas or its luminosity will scale with the dynamical mass of the accreting 
object. 

The observed upper limits reported here, therefore, show that the extreme object VIR- 
GOHI 21 is not an extreme X-ray object (in the sense of a high fx/ fo) nor has it accreted 
very much gas from its surroundings during its lifetime. 
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Fig. 1. — Swift XRT images of VI RGOHI 21, smoothed with a Gaussian kernel of <t=7" . Grey 
contours are reproduced from the iMinchin et al.l (120071 ) Hi observations with the Westerbork 
telescope. The black box represent the location of the putative Hi rotating disc. The circle 
centered at RA=12 7l 17 m 23 s , <5=14°39'29" represents the region from which we extracted the 
counts of the serendipitous source SWIFT-XRT J121723+1439.5. 
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Fig. 2. — (a): Upper limits to the mass of diffuse gas in VIRGOHI 21, estimated assuming 
a detection threshold of 17 counts in these XRT observations. The Swift XRT bandpass is 
not sensitive to emission from diffuse gas at temperature kT < 0.1 keV. Solid line : uniform 
density model, using a detection threshold of 8.2 counts in a 3.5' x 1.375' box. Dashed lines 
: (3 model with r c =0.5', using a detection threshold of 10.7 counts in a circle of radius 
1.75' (bottom line: /3=0.67; top line: /9=0.33). Given that the gas mass is proportional to 
the square root of the number of detected counts, the estimates are not very sensitive to 
changes in the detection threshold, (b) Upper limits to the luminosity of point sources in 
the VIRGOHI 21 field, estimated assuming a detection threshold of 4 source counts. The 
luminosity is proportional to the number of detected counts. 



